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Abstract. With the emergence of swept-volume ultrasound (US) probes, pre- 
cise and almost real-time US volume imaging has become available. This offers 
many new opportunities for computer guided diagnosis and therapy, 3-D images 
containing significantly more information than 2-D slices. However, computer 
guidance often requires knowledge about the exact position of US voxels relative 
to a tracking reference, which can only be achieved through probe calibration. 
In this paper we present a 3-D US probe calibration system based on a mem- 
brane phantom. The calibration matrix is retrieved by detection of a membrane 
plane in a dozen of US acquisitions of the phantom. Plane detection is robustly 
performed with the 2-D Hough transformation. The feature extraction process is 
fully automated, calibration requires about 20 minutes and the calibration system 
can be used in a clinical context. The precision of the system was evaluated to 
a root mean square (RMS) distance error of 1.15mm and to an RMS angular er- 
ror of 0.61°. The point reconstruction accuracy was evaluated to 0.9mm and the 
angular reconstruction accuracy to 1.79°. 



1 Introduction 

Until recently, 3-D ultrasound (US) volumes had to be manually reconstructed from a 
number of 2-D US slices acquired while slowly moving a 2-D probe over the target 
region. The so-called 3D freehand method is time-consuming, imprecise and not usable 
for many clinical applications requiring real-time acquisition [1|. The emergence of 
3-D swept-volume US probes solved most of the enumerated problems: a mechanical 
device capable of sweeping the 2-D crystal array of the probe over a target region makes 
it possible to acquire 3-D US volumes accurately and almost in real-time (Is to 4s per 
acquisition^. 

These new capabilities open an entire new field of applications in the domain of 
computer guided medical interventions based on US imaging. One can imagine tool 
guidance systems that would operate with permanently updated US volumes, visual- 
izing for instance slices at the tool tip position. More sophisticated applications could 



3 However, most currently available systems don't yet provide a real-time data transfer interface 
for 3-D data. Nevertheless one can acquire so called "4-D" images of three orthogonal volume 
slices in real-time using a video-capture device. In the rest of this article we make abstraction 
of this restriction, hoping that it will disappear with the next generation of 3-D echographs. 



carry out target localization inside the volumes through real-time registration and seg- 
mentation techniques, thus allowing to match pre-operative planning with intra-operative 
data. 

However, US-based guidance often requires knowledge about the position and ori- 
entation of the US volume in space. When using a tracking system this can be achieved 
by calibrating the US acquisition volume with a tracking reference fixed on the probe. 
Unfortunately it is virtually impossible to derive the calibration parameters directly 
from the geometry and parameterization of the probe. Almost all existing calibration 
systems rely therefore on statistical or segmentation-based object matching methods. 

1.1 Calibration Methods Overview 

A variety of techniques for 2-D US calibration was proposed in the literature; a com- 
prehensive review being given in 0. Calibration methods can be classified with respect 
to the target geometry they rely on. Single-point target methods identify a point, i.e. a 
bead, a calibrated pointer tip or a cross-wire, in the US image II3I4I5I . The difficulties 
consist in automatic geometry extraction in the US slice and US beam alignment with 
the phantom. 

Multi-point target phantoms are extensions of the single-point bead or cross-wire 
phantoms. They consist of a number of point targets with precisely known coordinates 
in phantom space. Their geometric configuration makes it possible to derive the calibra- 
tion parameters from the distances between the reconstructed intersection points visible 
in the 2-D US scan [6 |. Compared to single-point phantoms they require less image ac- 
quisitions due to their more discriminative geometry, but share the phantom alignment 
and feature extraction problems. 

Z-fiducial or N-fiducial phantoms address the alignment problem of point target 
methods. A calibration point is determined from the intersection points of a number 
of nylon strings with the US beam. This is possible due to a sufficiently discrimina- 
tive wire geometry [7 8 9|. Fiducial methods are more robust than point target methods 
but the difficulties concerning fully automatic feature extraction subsist. Also, Z- or 
N-fiducial phantoms require a high manufacturing accuracy to achieve a satisfying cal- 
ibration quality. 

Wall phantom methods are based on detection of the intersecting line of a planar 
surface with the 2-D US beam. In [4], a water tank bottom is imaged for calibration. 
The authors of [ 10] address the reverberation and line thickness problems inherent of 
wall phantoms by using a membrane variant. Both phantoms have difficulties when 
confronted with steep angles between the US beam and the plane because they cause 
line intensity and line sharpness degradation [2]. The Cambridge phantom scans a ro- 
tating bar, thus creating a virtual plane, to solve these problems |4|. The advantage of 
plane phantoms lies in the robustness of the feature extraction process which can, as a 
consequence, be reliably and fully automated. The pitfall of this method lies in the non- 
discriminative phantom geometry which can result in underdetermined systems if the 
acquired calibration samples do not cover all degrees of freedom. This can be avoided 
by strictly respecting the acquisition protocols presented in 1411 11 . 

Registration Phantoms: the last class of calibration methods relies on surface or 
intensity based registration techniques and therefore has the advantage of being inde- 



pendent of phantom geometry. The only requirement on phantom shape is that its US 
image is sufficiently discriminative with respect to rotations and translations, which 
is true for non-symmetric phantoms. The lack of precision of registration algorithms 
is the major drawback of this approach. To our knowledge, only one study examined 
registration-based 2-D probe calibration, registering US slices with an MRI image of 
the phantom 11121 . A 3-D approach is discussed in the next paragraph. 

1.2 3-D Probe Calibration 

Until today, only few studies about calibration of 3-D probes were carried out. Poon 
and Rohling lfl3l compared 3-D calibration based on a IXI-fiducial wire phantom, a 
pointer tip phantom and a cube phantom. The IXI wire phantom and the cube phan- 
tom methods require only one volume acquisition for calibration. The presented feature 
detection is semi-automatic. The best results yielded the IXI phantom with a mean er- 
ror in reproducibility of 1.5 mm, a RMS error of the point accuracy measure of 2.15 
mm and a RMS error of the reconstruction accuracy by distance measure of 1.52 mm. 
Bouchet et al lfT4l examined Z-fiducial phantom and achieved a RMS point accuracy 
error of 1.1mm. Two variants of a surface registration based 3-D calibration method 
were presented by Lange and Eulenstein in 1 15 1. The first one registers 3-D US images 
of the phantom with a geometric model derived from its CT scan. The second variant 
registers a number of US images of the phantom acquired from different positions. In 
both cases, surfaces are extracted manually. The authors claim that the latter approach 
could be fully automated. The CT variant performed best and yielded a RMS error in 
reproducibility precision of 1.8 mm and a RMS error in point accuracy of 2.0 mm. The 
ultrasound speed distortion problem is not addressed. 

In this study we propose a 3-D US calibration method based on a single plane mem- 
brane phantom. A fast, precise and accurate 3-D feature extraction algorithm relying on 
the 2-D Hough transform is presented. In contrast to existing 3-D US calibration sys- 
tems, the feature extraction process is fully automated. In the result section, precision 
and accuracy assessments are carried out using a specially designed validation phantom. 

2 Materials and Methods 
2.1 Acquisition Hardware 

The acquisition hardware consists of a GE Voluson 730 Pro 3-D US scanner and a 
NDI Polaris optical tracker with a 0.25 mm RMS error (as communicated by NDI). 
The tracking system operates with wireless (passive) infrared-reflecting rigid bodies 
equipped with flat markers. The ultrasound volumes are acquired with a 5 to 9 Mhz two 
dimensional curved array probe (see Fig.QJ). The piezo array of the probe is mounted 
on a mechanical device which is capable of sweeping regularly around its rotation axis 
within a predefined angular range. During the continuous sweeping process the US 
hardware reconstructs 3-D volumes from the series of acquired 2-D slices. The 3-D 
acquisition time ranges from Is to 4s, depending mainly on sweep angle and axial 
acquisition depth. Images are digitally transferred using a proprietary software from 
GE Medical Systems named 4D View. The US scanner also communicates the voxel 
size. The scan converter assumes the speed of sound (SoS) in tissue to be 1540 m/s. 



2.2 The Membrane Phantom 



The calibration phantom being dedicated to a clinical context, ergonomics considera- 
tions had an important impact on its design. We use a variation of the wall phantom 
presented in [4|, which is based on imaging the bottom wall of a water tank. The geo- 
metric form of the wall, which is a line in 2-D and a plane in 3-D, can be very robustly 
extracted from the US data using statistical algorithms like the Hough transform. This 
makes it possible to fully automate the feature extraction process without significant 
loss in precision and accuracy. This represents a big advantage over semi-automatic 
point-detection based phantoms in terms of calibration speed and ease of use. To over- 
come the plane thickness and the reverberation problems observable in US images of 
rigid surfaces [4| a filigrane nylon mesh membrane, tightly spanned on a planar rigid 
support with a circular and about 20cm wide hole, is used as target (see Fig.QJ). Rever- 
beration is further reduced by inclining the membrane plane with respect to the water 
tank bottom by 45°. A tracking reference (rigid body) is mounted on the membrane 
frame for phantom localization. The phantom is filled with water and equipped with a 
thermometer to measure water temperature. 




(a) (b) 



Fig. 1. Calibration hardware. Figure (a) shows an endorectal US probe mounted on 
an articulated arm. Figure (b) shows the membrane phantom. Both the probe and the 
membrane are equipped with infra-red reflecting passive rigid bodies for tracking. 



2.3 3-D Calibration Mathematics 

As illustrated in Fig. [2] four references are relevant for calibration: first of all, the mem- 
brane space M is defined as a reference in which the membrane lies in the origin and is 
parallel to e x and e y base vectors of M. In this space, every point with a zero z-ordinate 
is a membrane plane point. The phantom space Ph and the probe space Pr are de- 
fined by the rigid bodies that are attached on the phantom and on the probe. Finally, the 



US volume space U corresponds to the voxel space of the 3-D images acquired by the 
ultrasound device. Tph2M, Tp r 2Ph and Tu2Pr are homogenous 4x4 transformation 
matrices. 

Suppose that we identified a point p = (a;, y, z) in a US volume U as a point 



belonging to the membrane. With s 
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where Tph2M is known from membrane pre-calibration (see chap. \2A\ and Tp r 2Ph 
is given by the tracking system. Further, the scaling vector s is communicated by the 
US hardware. The remaining unknown element is the homogenous rigid transformation 
Tu2P- For convenience we define the elements of Tp r 2Ph ■ Tph2M as and the 
elements of Tu2P as bij G 1..4). The zero component of ([l} yields then 



= a 31 (s x xbn + s y yb\i + s z zbi 3 + b u ) + 
a32 (s x xb 2 i + s y yb22 + s z zb 2 3 + hi) + 
a33 (s x xb 31 + s y yb 3 2 + s z zb 33 + b 3i ) + 
a 34 . 



(2) 



Using Euler angles and a three-dimensional vector we can represent Tu2P with six 
variables, which leaves us in total with 6 unknowns to solve for. A detected plane can 
be added to the equation system by adding at least three plane points (Using of course 
the Tp r 2Ph measured while acquiring the US volume in which the plane was detected). 




Fig. 2. Illustration of the transformations involved in the calibration process. Note that 
scaling is omitted from the scheme for simplification. 



2.4 Membrane Pre-calibration 



To reduce the number of degrees of freedom of the calibration process, the membrane 
space M is determined using a pointer equipped with a tracking reference. A large 
number of surface points of the membrane-supporting structure is acquired in order to 
compute the plane equation using a least square approximation combined with a simple 
M-estimator to increase robustness. Since the phantom rigid body is permanently fixed 
on the phantom, pre-calibration has to be carried out only once. 

2.5 Acquisition Protocol 

The major drawback of single-plane phantoms resides in their barely discriminative ge- 
ometry. A plane can be described with only three variables, from which it follows that 
even with an optimal acquisition protocol, a minimum number of two acquisitions is 
necessary to cover all degrees of freedom. To obtain robust results the twelve-step ac- 
quisition protocol presented in (4) is used. The protocol improvement presented in ITD 
mainly addresses the z-axis imprecision problem inherent of most 2-D calibration sys- 
tems. Since 3-D probes give as much information in z-direction as in x- or y- direction 
this modification yields no particular advantage in the 3-D domain but requires at least 
18 steps. For that reason we stick to the original version. 

Sweeping and volume reconstruction being a continuous process of 1 to 4 seconds, 
significant distortions can be observed in the US volume when the probe is moved 
rapidly. Also, no direct access to the digital data is available which prevented synchro- 
nization of probe position measurement with US image acquisition. Therefore an artic- 
ulated arm for complete probe immobilization during acquisition is used, eliminating 
all motion-induced artifacts and time lags. Furthermore, immobilizing the probe makes 
it possible to perform high precision position measuring based on a large number of 
measures and outlier elimination. 

2.6 Feature Extraction 

The first step of the feature extraction process consists in correcting the distortion 
caused by the difference between US speed in water at room temperature and in human 
tissue at 37°. To determine US speed in water in function of temperature the polyno- 
mial formula established by Bilaniuk and Wong was used II16I171 . A distortion geome- 
try overview for all common probe types is given in [ 18 1. The distortion geometry of a 
sectorial probe is given in Fig.|3^. With v f w being the US velocity in water for a given 
temperature t and vt being the velocity in tissue, dx is determined using the following 
formula: 

dr = — • dw (3) 

v w 

Sectorial probe speed correction requires manual definition of the US origin and the 
scan head surface radius. A graphical user interface was developed for this purpose 
(see Fig.03). Origin and surface radius have to be defined only once during calibration. 
Plane detection can be carried out using the 3-D Hough transformation, but it would 
take several minutes to compute the result. Fortunately it is possible to determine the 
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Fig. 3. US speed correction. Fig. (a) illustrates the correction geometry of a sectorial 
probe. Ous is the probe origin, r the probe surface radius, dyy the point's P distorted 
distance from the probe surface in water, dr the corrected distance and P' the corrected 
point. Fig. (b) shows the probe mask used to determine US origin and probe surface 
radius. 



plane with good precision by simply extracting its intersection with two arbitrary vol- 
ume slices, using the 2-D Hough transform. To facilitate and to accelerate US speed 
correction the xy and zy planes passing through the scan head origin were used. The 
Hough transform implementation uses intensity accumulation and the following thresh- 
old sh for an image I: 

sh = max{i <E Hist(I)} + (max{i el} - min{i E I})/3 . (4) 

The purpose of sh is to ignore the low-intensity water background, which represents 
the largest part of the image. 




(a) (b) 



Fig. 4. Screenshot of a successful automatic plane extraction. Note that lines are cor- 
rectly detected in spite of a degraded membrane image caused by a steep scan angle. 
The arrows point at the line intersections with the mask. 



2.7 Optimization 



Optimization of (f2]i is carried out with the non-linear Levenberg-Marquardt implemen- 
tation given in [ 19 1. A random restart scheme within a range of reasonable initialization 
values robustifies this process. 

2.8 Visual Back-tests 

The plane coordinates as resulting from the optimization process are visualized as a line 
in the slices used for feature extraction, and distance plus rotational errors between the 
segmented and the calculated line are evaluated. This allows for manual replacement of 
evident outliers with new acquisitions and for recomputation of the calibration without 
requiring a complete restart. 

3 Experimental Results 

3.1 Test Configuration 

Precision and accuracy assessments were carried out using a membrane plane pre- 
calibration with an RMS surface distance error of 0.43mm for the measured surface 
points. A total of ten calibrations were performed using the twelve-step protocol. The 
probe rigid body was not moved between calibrations. The probe was mounted on an 
articulated arm and immobilized during position and image acquisition. The water tem- 
perature was 23°. The acquired US volumes had a size of (199, 199, 199) isotropic 
voxels with 0.477mm side lengths. 

3.2 Feature Extraction Quality 

The Hough transform extracted lines correctly for 238 out of 240 acquired images. 
Line detection failed for images on which only a very small part of the membrane 
was visible. In these cases, lines had to be manually determined. Note that no manual 
outlier elimination was carried out. As detection failures were rare, calibration was in 
average carried out in about twenty minutes. To get a better idea about the quality of 
feature extraction and about the presence of distortions in the membrane images we 
measured the detection precision: using (O we can calculate the distance error between 
a measured plane point and the computed plane as follows: 

e(x, y, z) = a 31 (s x xbn + s y yb 12 + s z zb 13 + bu) + 

032 (s x xb 2 i + s y yb 2 2 + s z zb 23 + b 2i ) + 

033 (s x xb 31 + s y yb 32 + s z zb 33 + b 34 ) + 
a 34 . 

For each calibration, the average and the root mean square (RMS) distance of a set of 
points to the pre-calibration plane was computed using (0. For each line we computed 
ten equidistant points between the extreme points on the line segment inside the US vol- 
ume. The angular feature extraction error is defined as the angle between the computed 



plane normal and the cross product of the directional vectors of the two extracted lines. 
Based on this definition the maximum and the RMS angular errors were computed for 
each acquired volume of the calibration. The aggregated errors for all calibrations can 
be found in Table Q] 

Table 1. Aggregated Feature Extraction Precision. 



Distance Error Distance Error Angular Error 

[mm] [vox] [deg] 

RMS Error 0.37 0.77 0.26 

Max Error 1.30 2.73 1.09 



3.3 Calibration Precision 

The calibration precision measures the reproducibility of calibration results. Again, 
both the translational and the angular errors were assessed. The translational error is 
defined as the standard deviation of the volume center after scaling and right-hand 
multiplication to the different calibration transformations Tjj 2Pr . The angular error is 
measured as the standard deviation of angular differences between the (0, 0, 1) vector 
after scaling and right-hand multiplication to the different calibration transformations 
T|j 2 p r (see Tabled. 

Table 2. Calibration Precision. 



Distance Error Distance Error Angular Error 

[mm] [vox] [deg] 

Standard Deviation 1.15 2.41 0.61 

Max Error 1.99 4.03 1.12 



3.4 Reconstruction Accuracy 

Reconstruction accuracy was assessed using the bead phantom illustrated in Fig. [5] 
Note that the beads are co-planar within a precision of 0.25mm (RMS). The left-hand 
three beads form the left triangle while the right-hand beads form the right triangle. 
The distance du of the triangle bary centers was evaluated with an estimated accuracy 
of about 0.5mm. 

Twenty US images of the phantom were acquired, ten imaging the left triangle and 
ten the right triangle. Note that images of one triangle did not intersect with the sec- 
ond one. The bead centers were manually extracted from the images. The reconstructed 
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Fig. 5. Reconstruction Accuracy Measurement Phantom. 



triangle barycenters and normals were then projected into probe space for each calibra- 
tion, which yielded 100 point and vector pairs. The distance error for a point pair is 
defined as the difference between their Euclidean distance and ds- The angular error 
for a vector pair is defined as the angle between both vectors. The results are given in 
Table 



Table 3. Reconstruction Accuracy Results. 





Distance Error 


Distance Error 


Angular Error 




[mm] 


[vox] 


[deg] 


RMS Error 


0.90 


1.90 


1.79 


Max Error 


2.44 


5.11 


3.03 



4 Discussion 

Probe calibration for currently available swept-volume 3-D probes makes only sense for 
static applications. Real-time access to acquired 3D images is currently not provided by 
the Voluson hardware. Also, depending on scan parameters, the duration of the volume 
sweeping process ranges from Is to 4s. Due to probe or tissue motion, the physical 
location of voxels can therefore be way off the position indicated by the calibration. The 
latter problem could be reduced by scan-heads equipped with high-frequency sweeping- 
devices, but it will disappear only with non-sweeping 2-D piezo array US probes. Until 
then, it is better to calibrate probes using an articulated arm for immobilization. 

Passing from 2-D to 3-D calibration improves calibration results (for information 
on 2-D precision of the membrane phantom see [20|) because the z-axis-uncertainty 
inherent of 2-D calibration is eliminated: when giving a plane instead of a line as input 
to the optimizer the rotational degrees of freedom are significantly better covered. This 
allowed us to reduce the number of acquisitions for calibration while still achieving 
precise and accurate results. 



Feature extraction from membrane phantom images showed both robust and precise 
results. On our set of test images the line blurring and intensity degradation effects oc- 
curring when scanning a wall phantom from an oblique angle were correctly handled 
by the Hough transform: lines were consistently placed in the center of the beam width. 
Note that the feature extraction precision RMS and maximum errors reported in Ta- 
ble Q] are relatively small, which indicates that the physical plane location corresponds 
indeed with the beam width center line. Membrane reverberation was not observable 
and did therefore not disturb the detection process. Also, the membrane phantom was 
not exposed to line thickness problems. Due to these characteristics, feature extraction 
could be fully automated (up to the manual US origin and probe radius determination 
required for US speed correction). 

The user-independency resulting from automated feature extraction is partially coun- 
ter-balanced by the necessity to follow a protocol for data acquisition, which can re- 
introduce user bias. Nevertheless we believe that it is more convenient to follow a 
simple acquisition protocol instead of extracting features semi-automatically from US 
volumes. Further we preferred to correct for US speed errors instead of requiring 50° 
water or phantom fill materials that have the same US characteristics than human tissue, 
which makes it possible to use the phantom in a clinical context. 

The overall calibration time of about 20 minutes is mostly due to manipulation of 
the articulated arm, high precision probe position measurement and data transfer, which 
requires several manual interventions. Feature extraction and optimization is computed 
in several seconds. In cases where the feature extraction precision evaluation shows 
poor results, a visual verification and eventually a correction have to be carried out, 
which requires some additional minutes. Significant speed-up could be achieved by 
automating the communication between the US scanner and the calibration computer. 

The presented calibration system assumes that the SoS in human tissue is uniform 
and that it corresponds to the SoS internally used by the US scanner, which is in general 
the mean SoS in tissue of 1540m/s. However, SoS varies with different types of tissue: 
The SoS in fat is approximately 1450m/s, in blood 1570m/s, in the brain 1541m/s and 
in water 1480m/s. As the US generally crosses tissue layers of different thickness and 
different types on its way through the body, and as the target tissue is often viewed from 
various positions, the in vivo accuracy of the calibration may show fluctuations of more 
than 5 per cents in extreme cases during an examination. Also, for some applications it 
would be appropriate to use a different mean SoS than the 1540m/s for calibration, but 
this is beyond the scope of this study. 

Future work will address the twelve-step acquisition protocol which contains a lot of 
redundancy. Also, our system does currently not provide a foolproof indicator for miss- 
ing coverage of degrees of freedom. We therefore started experiments with an Eigen- 
value system similar to the one presented in |2D . 

5 Conclusion 

A robust 3-D US probe calibration method designed for clinical usage was presented. 
Calibration can be carried out in about twenty minutes due to fully automatic 3-D Plane 
extraction based on robust and efficient 2-D line detection. The point reconstruction 



accuracy of our phantom can compete with previously presented 3-D phantoms: Lange 
and Eulenstein communicated RMS errors between 2.0mm and 2.2mm 1 15], Bouchet 
et al were confronted to 1 .1mm RMS point accuracy [ 14] while Pohn and Rohling pub- 
lished errors between 1.52mm for their IXI-wire, 1.59mm for the cube and 1.85mm 
for their stylus approach. With 0.90mm RMS point accuracy (see TableO we achieved 
slightly better results. Finally, the proposed method is temperature-independent and 
uses water as transmission matter which facilitates its usage. 
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